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We propose a chiral symmetry restoration mechanism in monolayer graphene, in analogy with the 
strongly coupled gauge theory. The chiral (sublattice) symmetry of graphene, which is spontaneously 
broken under the effectively strong Coulomb interaction, is restored by introducing the Kekule- 
patterned lattice distortion externally. Such a phase transition is investigated analytically using 
the techniques of strong coupling expansion on the lattice gauge theory model, by preserving the 
original honeycomb lattice structure. We discuss the relation between the chiral phase transition 
and the spectral gap amplitude, and we show the modification of the dispersion relation of the 
quasiparticles. 



PACS numbers: 73.22.Pr,71.35.-y,11.15.Ha,11.15.Me 

Ever since its first experimental isolation in 2004 1 , 
graphene has attracted a great deal of interest in the 
fields of particle physics as well as materials science^. 
One of its important features is that the electrons/holes 
on monolayer graphene can be effectively described as 
massless Dirac fermions around half filling 3,4 . Since 
the fermions on the (2+l)-dimensional sheet interact 
with the electromagnetic field in the (3+l)-dimensional 
space, one can apply "reduced quantum electrodynamics 
(QED)"— as an effective field theory for this real material. 

In reduced QED, the chiral symmetry is suggested to 
be spontaneously broken if the Coulomb interaction be- 
tween fermions is effectively strong. In graphene, such a 
strong Coulomb interaction is supposed to be achieved 
if the electromagnetic field is not screened by dielectric 
substrates, i.e. the layer is suspended in the vacuum. 
Breaking of the chiral symmetry corresponds to the spon- 
taneous breaking of the inversion symmetry between two 
triangular sublattices, leading to the generation of a mass 
gap£. This mechanism is analogous to the dynamical 
mass generation mechanism of quarks in quantum chro- 
modynamics (QCD)-. The techniques of lattice gauge 
theory have been employed, on the hypothetical square 
lattice, to investigate the chiral symmetry breaking in 
graphene either numerically^ or analytically^. 

In this paper, we propose a mechanism that restores 
the "chiral symmetry" even under the strong Coulomb 
interaction. Such a phase transition is induced by the 
external Kekule distortion, which is described by the 
alternating bond strength like the benzene molecule 1 ^. 
The Kekule distortion yields the asymmetry between two 
Dirac points and gives a finite gap in the spectrum, but 
does not break the chiral (sublattice) symmetry explic- 
itly. It can be introduced externally by the effect of the 
substrate or the addition of atoms on the laye r 11 ' 12 . The 
interplay between the breaking of the chiral symmetry 
and the Kekule distortion at long wavelength, i.e. in 
the mean field, has not yet been understood in previous 
works 13 . 

In order to incorporate the effect of the Kekule dis- 
tortion, the effective gauge theory model with the orig- 
inal honeycomb lattice structure would be helpfu l 14 ' 15 . 
First we construct a U(l) lattice gauge theory model for 



graphene, with both the spatial and temporal directions 
regularized on the hexagonal-prismatic lattice. Then, we 
apply to this model the strong coupling expansion tech- 
niques of the lattice gauge theory, which has been em- 
ployed widely in analyzing the properties of strong cou- 
pling gauge theories such as QC D 16 ' 17 . We demonstrate 
the chiral symmetry restoration transition by adding the 
Kekule distortion term to the model action. The inter- 
play between two orders results in the modification of the 
gap size, depending on the distortion amplitude. 

Noninteracting electrons on the monolayer graphene 
with the perfect honeycomb lattice symmetry are de- 
scribed by the conventional tight-binding Hamiltonian 



H = h H E 

rg J 4i=l,2,3 



a t (r)6(r + s i ) + & t (r + s I )a(r)] , (1) 



with the nearest-neighbor hopping amplitude h. The op- 
erators (a) and (b) create (annihilate) a fermion on 
the triangular sublattices A and B, respectively. Sj (i = 
1,2,3) is the orientation vector pointing from an A-site 
to a neighboring B-site, with the lattice spacing |sj| = 
a Hc = 1.42 A. This Hamiltonian yields the energy eigen- 
value E(k) = ±/i|$(k)|, where $(k) = £i=i,2,3 e*' s «. 
E(k) vanishes at two Dirac points K± in the Brillouin 
zone SI (see FigQJ. By neglecting quadratic and higher 
order terms in momentum, $ can be written around the 
Dirac points as $(K±+k) = (3/2)a Hc (±k x +ik y )+0(k 2 ), 
which we call here the "linear approximation" . Thus, the 
dispersion relation can be linearized around the Dirac 
points, with the Fermi velocity v F = (3/2)a Hc /i about 
300 times smaller than the speed of light 3 . An effective 
field theory description has been established in contin- 
uum, with a 4-component Dirac fermion ip = [a(K + + 
k), a(K_ + k), b{K+ + k), 6(K_ + k)] T 4 

From the tight-binding Hamiltonian in Eq.([T]), we can 
construct a Euclidean action for fermions with the orig- 
inal honeycomb lattice structure. As done in Ref0, we 
perform a temporal scale transformation r — > t'/v f , so 
that the effective Fermi velocity in the action can be ren- 
dered to be unity. Here we discretize the temporal di- 
rection with the lattice spacing a T i, which is set to be 
comparable to the spatial lattice spacing a Hc . (In this 
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by ft, which is 0.037 in the vacuum-suspended graphene. 
(Here we take v F ~ c/300 at any value of (3.) In 
this work, we only consider the leading order [0(/3°)] in 
the strong coupling expansion, which corresponds to the 
strong coupling limit ((3 — s> 0) of the Coulomb interaction. 
The gauge term Sg does not contribute in this limit, so 
that the partition function of this system can be written 
only in terms of the fermionic term Sf as 



FIG. 1: A schematic picture of the Brillouin zone correspond- 
ing to the honeycomb lattice. The Brillouin zone Q (gray 
rhombic region) is spanned by the reciprocal lattice vectors 
Ki and K2. If the Kekule distortion pattern [Eq.([7|] is in- 
troduced, Q is split into three hexagonal cells: Cl and Cl±, 
surrounding k = and K± (Dirac points) respectively. 



paper, we set a T i equal to a Hc for convenience of calcula- 
tion.) A pair of fermion doublers are inevitably generated 
in the temporal directioniS,, which we reinterpret as the 
spin (up/down) degrees of freedom. 

In order to incorporate the interaction between the 
fermions and the electromagnetic field, we require the 
local U(l) gauge invariance, with the compact link vari- 



[daUa][dbUb][dU T ,]e 



(3) 



exp 



te J^ +a ^' dr'Ai and £7;(r, r') 



ables U T /(r,T / ) 
exp ~ 

obtain a "gauged" "honeycomb lattice model 
1 



f* s ' dr ■ A (e is the electric charge); thus we 
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- [a) (x)U T '(x)a(x + a T >) 

reA;r' 
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[b\x)U T '(x)b{x + a T > 



H.c.] 
H.c] 



(2) 



First we observe the behavior of the global symme- 
try of the system when there is no lattice distortion. 
By the integration over the link variables U T >, tem- 
poral kinetic terms [first two lines in Eq.@] are con- 
verted into 4-Fermi terms — (1 / 4)n a (x)n a (x + a T <) and 
— (1 / 4)nb(x)ni,(x + a T i), which are spatially local and 
temporally non-local by one lattice spacing. Here the 
bosonic operator n x (x) = x^{ x )xi x ) (x — a ,b) corre- 
sponds to the local charge density. These terms can be 
considered as the on-site (Hubbard) interaction, in which 
fermions on the same position with opposite spins inter- 
act with each other. Such an interaction may lead to the 
spontaneous breaking of the sublattice symmetry at the 
tree level, as long as the spin symmetry is not broken. 

Here we introduce an auxiliary field a = (n a — fib) 
by the Stratonovich-Hubbard transformation, which cor- 
responds to the charge density imbalance between two 
sublattices. It corresponds to the "chiral condensate" 
(ipip), the order parameter of chiral (sublattice) symme- 
try breaking, in the continuum effective theory. By the 
mean-field approximation over a, the effective action at 
zero temperature is given in the quadratic form 



+ - 



a T >h 



r£A,r' i=l 



where x denotes the (2+l)-dimensional position (r,r'). 
Here we take a) and W (a and b) as the Grassmann fields 
corresponding to the creation (annihilation) operators 
of the fermionic quasiparticles around half-filling. Ki- 
netic term of the gauge field, Sg, can also be discretized 
on the honeycomb lattice with the U(l) link variables. 
Sq becomes proportional to the inverse of the effective 



, . (4) 

where 'F = (a,b) T and a = a T >h/v F . Thus the inte- 
gration over the fermions can be successfully performed, 
yielding the free energy per one pair of A- and B-sites in 
the strong coupling limit, 



Coulomb coupling strength, 



r /e 2 = 4we/c 



QED ' 



where e is the dielectric constant of the surrounding ma- 
terial and a QED ~ 1/137 is the fine-structure constant of 
QED. 

This lattice model reproduces the "reduced QED" -like 
model 19,20 in the continuum limit. The fermions (quasi- 
particles) propagate much slower than the electromag- 
netic field (photons) (i.e. v F <C c), so that the retarda- 
tion (magnetic) effect of the electric field can safely be 
neglected. Therefore, the spatial link variables Ui can 
be set to unity here, which we call the "instantaneous 
approximation" . 

Now we can perform the strong coupling expansion 



(a) is dominated by the second term (singularity 
from fermion one-loop) around a = 0, while it is domi- 
nated by the first term (tree level of the auxiliary field) 

at large |er|. Therefore, the effective potential F^(a) 
has a minimum at finite | cr] = 0.343, which gives the 
charge density imbalance between two sublattices. Finite 
a breaks the sublattice symmetry of the honeycomb lat- 
tice, leading to a gap opening in the fermionic spectrum. 
The term (o /2)(a) a— b^b) in Eq.Q, which corresponds to 
the effective mass term m^ipip in the 4-component Dirac 
fermion description, modifies the dispersion relation into 
E(k) = ±V|4>(k)| 2 + (<; F (7/2a r 2 . 
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|a$(k) 
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FIG. 2: The free energy per one pair of A- and B-sites in the 
strong coupling limit, Fjg , as a function of a. "Honeycomb" 
shows -Feff m Eq.([5]l with the exact form of $(1^). "Lin- 
ear" shows F^ff (a) in Eq.© with the linear approximation 
$(K±+k) ~ (3/2)a Hc (±k x + ik y ). "Square" shows the effec- 
tive potential obtained by the strong coupling expansion with 
the square lattice formulation^ per four square lattice sites, 
corresponding to a pair of honeycomb lattice sites. 



In FigJU we compare the result from the honeycomb 
lattice formulation with that from the square lattice for- 
mulation given in Ref§. Here, F?' from the square 
lattice formulation is given by the free energy per four 
square lattice sites, which corresponds to a pair of A- and 
B-sites on the honeycomb lattice. We also show F^\a) 
in Eq.© with the linear approximation of ^(k) in Figj2] 
Since they have the Dirac cone structure around k = K± 
in common, the logarithmic divergence around a = ap- 
pears independently of the lattice formulation, leading 
to the spontaneous sublattice symmetry breaking. The 
quantitative difference in F^g comes from the difference 
of the dispersion relation in the quadratic and higher or- 
der terms in momentum, which are neglected in the linear 
approximation but become dominant far from the Dirac 
points in the momentum space. 

Let us now examine the effect of an external Kekule 
distortion on the sublattice symmetry; it is represented 
by the additional Hamiltonian 



Hk = J2J2 Mn(r)at(r)&(r + s;) + H.c. 



r£A i=l 



i(K + -s;+G-r) 



i(K_- 



G 



(6) 



(7) 



where G = K + — K_, and A denotes the amplitude 
of the external Kekule distortion 21 . In this paper, we 
consider the case in which the distortion amplitude A is 
real and uniform. Since the distortion pattern is peri- 
odic by three times the unit cell of the honeycomb lat- 
tice, the Brillouin zone fl is split into three hexagonal 
regions, fl and f2±, surrounding k — and k = K± 
respectively (see FigJT]). For free fermions, the Kekule 
distortion opens a finite gap in the dispersion relation of 
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FIG. 3: Charge density imbalance between two sublattices, a , 
and the total energy gap E(K±; A) = (v F a/2a T /) 2 + A 2 , 
as the functions of the external Kekule distortion A. <j(A) 
first grows and then decreases as |A| increases, vanishing at 
the critical value Ac/Sh = 0.485, while E(K±;A) remains 
finite for any value of A. The dotted line shows the gap 
amplitude of the free fermion. 



quasiparticles, E hcc (K± + k; A) = y^kp + |A| 2 + 0(fc 4 ), 
while it preserves the sublattice symmetry of the honey- 
comb lattice. The magnitude of the gap is proportional 
to the distortion amplitude A. In the Dirac fermion de- 
scription, Kekule distortion can be written as an exter- 
nal field A?/>73?/;, which is invariant under the U(l) chiral 
transformation generated by the chirality (sublattice in- 
version) 75. 

By adding the Kekule term to the lattice effective ac- 
tion in Eq.Q, it is modified as 

F ~ ^ 4 + ^ 1 ' \ a<f (k) -(a/2)lj ( ' >: 



reAUB-T' kGO:T' 



(8) 

where *(k) = [o(k), a(K+ + k), a(K_ + k), b(k), b(K+ + 
k), 6(K_ +k)] T and ^3 is a 3 x 3 unit matrix. <J> is defined 
by a 3 x 3 matrix, 

/ $(k) A$(K_ +k) A$(K+ +k) 
$(k)= A$(K_+k) $(K++k) A$(k) 

\A$(K + +k) A$(k) $(K_+k) 

(9) 

with the dimensionless distortion amplitude A = A/3h. 
By integrating out the fermionic degrees of freedom, we 
obtain the effective potential F^g as a function of er, 
with the external parameter A. When there is no Kekule 
distortion, it reproduces the effective potential in Eq. ([5]) . 

The minimum of the effective potential f|j (<t; A) yields 
the expectation value of a as a function of A, as shown 
in Figj3] By diagonalizing the matrix in Eq.© at k = 0, 
we obtain the magnitude of the spectral gap as 



E(K±; A) = y [v F o{A)/2a T >] + A 2 , 



(10) 



which is modified from the gap of the free electrons by 
the distortion-dependent exciton condensate <r(A). 

As seen from FigSl G grows around |A| = and de- 
creases at large |A|, with the critical value |Ac/3/i| = 
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0.485. The reduction of cr(A) can be qualitatively un- 
derstood in terms of the effective field theory with a 4- 
component Dirac fermion: the amplitude of the fermion 
propagator (k-'f+a/2+Aj 3 )~ 1 gets suppressed as | A| — > 
oo, leading to a reduction of the dominance of the fermion 
one-loop effect in the effective potential. The sublattice 
symmetry is fully restored, i.e. a vanishes, at the criti- 
cal value Ac, where <r(A) exhibits the mean- field critical 
exponent z = 1/2 with the definition a oc |A — Ac| z . 
As a result, we find that the external Kekule distortion 
can restore the sublattice symmetry of the honeycomb 
lattice, although the gap in the spectrum remains finite, 
as shown in Figj3] On the other hand, c(A) increases 
quadratically as cr(A) = 0.343 + 7.02|A| 2 + 0(A 4 ) for 
small |A|. It comes from the quadratic and higher order 
terms in momentum which are neglected in the contin- 
uum effective theory; if the linear approximation is ap- 
plied to F^\a; A), the resulting ct(A) decreases mono- 
tonically around A = 0. 

In this paper, we have investigated the interplay be- 
tween two symmetry breaking patterns of the honeycomb 
lattice, in the strong coupling limit of the Coulomb inter- 
action. Due to the on-site part of the Coulomb interac- 
tion, the sublattice symmetry is spontaneously broken 
if the system has the perfect honeycomb lattice sym- 
metry. As we introduce the Kekule distortion as an 
external field, the system reveals a second-order phase 
transition - restoration of the sublattice (chiral) sym- 
metry. As a result, the gap amplitude of monolayer 
graphene in the strong coupling limit is modified from 
that of free fermions if the external Kekule distortion 
is under the critical value, while it is not modified at 
large Kekule distortion. If the "chiral symmetry bro- 



ken" phase is achieved experimentally in the vacuum- 
suspended graphene, the modification of the gap ampli- 
tude will have an effect on gap engineering, which has 
recently become important for industrial applications as 
electronic devices. 

While the Kekule distortion is introduced as an 
external parameter in the present paper, it is proposed 
to occur spontaneously if the interaction between the 
nearest-neighboring (NN) or the second-NN sites is taken 
into account22r— . We are currently studying such an 
effect in the next-to leading order [C^/? 1 )] in the strong 
coupling expansion, which includes the NN interaction 
Z)j=i [a t (a;)&(x + Sj)tf(x + Sj + a T >)a(x + a r >) + H.c] . 
With a sufficiently strong NN interaction, our pre- 
liminary study shows a phase transition between the 
sublattice symmetry-broken phase and the spontaneous 
Kekule distortion phase. Such a phase transition is 
expected to be of first order, since either the chiral 
condensate or the Kekule distortion remains finite due 
to the logarithmic behavior in the effective potential. 
The effect of the phonon-mediated interaction is also of 
great importance; if such an interaction is taken into 
account, the superconducting order may also have to be 
considered^^. 
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